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ABSTRACT 



We report Hubble Space Telescope imaging, obtained 155 and 449 days after 
the 2006 outburst of the recurrent nova RS Ophiuchi, together with ground- 
based spectroscopic observations, obtained from the Observatorio Astronomico 
Nacional en San Pedro Martir, Baja Cahfornia, Mexico and at the Observato- 
rio Astrofisico Guillermo Haro, at Cananea, Sonora, Mexico. The observations 
at the first epoch were used as inputs to model the geometry and kinematic 
structure of the evolving RS Oph nebular remnant. We find that the modeled 
remnant comprises two distinct co-aligned bipolar components; a low-velocity, 
high-density innermost (hour glass) region and a more extended, high-velocity 
(dumbbell) structure. This overall structure is in agreement with that deduced 
from radio observations and optical interferometry at earlier epochs. We find that 
the asymmetry observed in the west lobe is an instrumental effect caused by the 
profile of the HST filter and hence demonstrate that this lobe is approaching the 
observer. We then conclude that the system has an inclination to the line of sight 
of 39;Jo de grees. This is in agreement with the inclination of the binary orbit and 
lends support to the proposal that this morphology is due to the interaction of 
the outburst ejecta with either an accretion disk around the central white dwarf 
and/or a pre-existing red giant wind that is significantly denser in the equato- 
rial regions of the binary than at the poles. The second epoch HST observation 
was also modeled. However, as no spectra were taken at this epoch, it is more 
difficult to constrain any model. Nevertheless, we demonstrate that between the 
two HST epochs the outer dumbbell structure seems to have expanded linearly. 
For the central (hour glass) region there may be evidence of deceleration, but it 
is harder to draw firm conclusions in this case. 
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Subject headings: linerprofiles — binaries: symbiotic — circumstellar matter — novae, 
cataclysmic variables — stars: individual (RS Ophiuchi) 
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INTRODUCTION 



RS Ophiuchi is a symbiotic recurrent nova (RN) with outbursts r ecorded in 1898, 19 33 



1958, 1967, 1985 (see 



Rosino 



1987 



Rosino Sz lijima 



19871 ) and 2006 flEvans et al. 



2008 



and references therein). However, during the period 1898 - 1933, befor e the star was 



known 



to be a RN, some outbursts might have occur red, most hkelv in 1907 



there is also evidence for an eruption in 1945 (jOppenheimer fc Mattei 



observed to be undergoing its latest eruption on 2006 February 12.83 (INarumi et al. 



Schaefer 



2004 ) and 



19931). RS Qph wa s 



20061 ). 



reaching a magnitude V = 4.5 at this time. We define this as t = 0. The opt ical light curv e 



then began a rapid decline, consistent with that seen in previous outbursts (IRosino 



0.3 kpc flBode 



1987, 



AAVS djJ). The dist ance to RS Qph has been derived from several lines of evidence as 1.6 ± 



1987: see also 



Barry et al. 



20081 and the discussion in Section 3.1 below). 



The RS Qph binary system compri ses a red giant, with a spectral type estimated 



to be in the range G5III to M4III (e.g. 



see 



Bohigas et al 



1989 : 



Anupama fc Mikolajewska 



1999), and a white dwarf 



(WD) 



Dobrzycka fc Kenyon 



1994 



with mass thought to be close to the Chandrasekhar lim it 



Shore et al 



1996 



Fekel et al 



2000 



Brandi et al 



20091 ). 



Accretion of hydrogen-rich material from the red giant onto the WD surface leads to a 
thermonuclear runaway, as in the outbursts of classical novae (CNe). The much shorter 
inter-outburst period for RNe compared to CNe is predicted to 
of high WD mass and a high accretion rate (e.g 



is predicted to 


De due to a combination 


Starrfield et al. 


1985; 


Yaron et al. 


2005) 



These models lead to the ejection of smaller amounts of material at higher velocities than 
those for CNe (typically 10~^ to 10^^ Mq and several thousand km s~^, respectively, for 
RNe). Spectroscopy of RS Qph has indeed shown Ha line emission with FWHM = 3930 



km s~^ and FWZI = 7540 km s~^ on 2006 February 14.2 {t = 1.37 days; 



Bull 



20061). 



^see http:/ /www.aavso.org 
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Unlike CNe, where the mass donor is a low mass main sequence star, the presence of 
the red giant in the RS Oph binary means that the high-velocity ejecta run into a dense 
circumstellar medium in the form of the red giant wind, setting up a shock system. The 
forward shock mov ing into the wind has a gas temperature ~ 2.2x10^ K for a shock velocity 



4000 km s"^ flBode et al. 



20061) 



Very Long Baseline Interferometry (VLBI) radio observations by 
at t = 13.8 days showed a partial ring of non-thermal radio emission resu 
expan ding shock, which developed into a more bipolar structure (see also 



O'Brien et al 



(120061 ) 



ting from the 



O'Brien et al. 



20081 ). They suggested that the asymmetry in this ring could be due to absorption in the 



overlying red giant w ind and they also note d the emergence of more extended components 



to the east and west. 



Sokoloski et al. 



( l2008l ) followed the evolution of the radio source with 
the Very Long Baseline Array (VLB A) at 4 epochs between 34 and 51 days after outburst. 
They found a central thermally-dominated source linked by what appeared to be a narrow 
(collimated) outflow to expanding non-thermal lobes which they interpreted as the working 
surfaces of jets. Jet colli mation could be du e to the expected accretion disc around the WD. 



It should be noted that 



Taylor et al. 



(119891 ) had also interpreted their 5 GHz VLBI map 
from day 77 after the 1985 outburst in terms of a central thermal source with expanding 
non-thermal lobes. 



From Hubbl e Syace 



( 12007I . hereafter 



Telescope (HST) observations, 155 days after outburst. 



Bode et al 



Paper j ) revealed the remnant to have a double ring structure with the 



major axis lying in the east- west direction (see also Section E]). They suggested that there 



was evidence for deceleration in the north- south direction by comparing t heir observations 

J 



to earlier epoch observations in the radio (lO'Brien et al. 



20061 ). 



Paper 1] also provided 



preliminary models of the remnant as a bipolar structure which implied a true expansion 
velocity, Vexp = 5600 ± 1100 km s~^ for the material at the poles. It was proposed that 
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the shaping of the remnant occurred due to the ejecta interaction with the pre-existing 
circumstellar environment. 

The projected 2-D geometry of nova explosions on the sky is usually complex, as are 
their 1-D spectra. However, taken together they provide the underlying 3-D structure of an 
object. The projected image on the sky provides spatial information on its structure while 
the spectra provide information about the component of the velocity vector along the line of 
sight. Here we re-examine the i/ST observations of the expanding nebular remnant taken at 
t = 155 days and report those at a second epoch at t = 449 days after outburst. We compare 
these to spectral observations (where available), then using a modeling code, determine the 
detailed 3-D structure of the remnant. In Section [2] we present our observations and data 
reduction methods. The results are presented in Section [3] and then modeled in Section HI 
In Section [5] we discuss the implications of the results found in Section |H Finally, in Section 
Owe summarise our conclusions. 

2. OBSERVATIONS AND DATA REDUCTION 

We imaged RS Oph with the HST on two occasions, under the Director's Discretionary 
programs GO/DD-11004 and GO/DD-11075. The first observations were made on 2006 
July 17 (155 days after outburst) with the High-Resolution Channel (HRC) of the Advanced 
Camera for Surveys (ACS), with a scale of 0'.'025 pixel"^ (equivalent to 40 AU pixel"^ at 
d = 1.6 kpc). The second observations were made on 2007 May 7 (449 days after outburst) 
with the Planetary Camera CCD of the Wide Field Planetary Camera 2 (WFPC2), with 
a scale of 0'.'046 pixel"^, since in the meantime the ACS had suffered a failure. During the 
first epoch observations, the ACS/HRC was used with three narrowband filters to isolate 
the Ha+[N ll] (F658N), [O ill] A5007 (F502N) and [Ne v] A3426 (F344N) nebular emission 
lines. For the subsequent WFPC2 observations, only the F502N filter was used in order to 



- 8 - 



isolate the [O ill] A5007 emission line. This line was chosen based on the results from the 
first epoch observations. 

Here we concentrate on the F502N filter observations common to both epochs. The 
other two HST filter bands, F344N and F658N, will not be discussed because we have 
no spectral information available for the former a nd the latter presented little significant 



structure in the epoch 1 HST images (see 



Paper j ) 



All data were reprocessed using standard procedures out 



WFPC2I3 Data Handbooks and the Pvdrizzlc 



were generated using TinyTim (IKrist 



ined in the AC^^ 



and 



and Multidrizzlcl Handbooks. PSF profiles 



19951 ) for both ACS/HRC and WFPC2 images. In 
the ACS/HRC and WFPC2 images, the PSF fiux was scaled to an extent such that when 
subtracted from the observed image the new PSF-subtracted image pixel fiux was close 
to those surrounding it. Deconvolution was then performed using the Lucy- Richardson 
method for the ACS/HRC and WFPC2 images. CLEAN and Maximum Entropy techniques 
were also used as a check, and produced similar results. 

Spectral observations were carried out at the Observatorio Astronomico Nacional en 
San Pedro Martir, Baja California, Mexico, using the 2.1 m telescope with the Boiler & 
Chivens and Echelle spectrographs, and at the Observatorio Astroffsico Guillermo Haro, 
at Cananea, Sonora, Mexico, using the 2.12 m telescope with their Boiler & Chivens 
spectrograph. The Boiler & Chivens spectrographs have intermediate spectral resolutions, 
from R ~ 1500 to 3500, while the Echelle spectrograph has a maximum resolution of about 



^See http://www.stsci.edu/hst/acs/documents/handbooks/DataHandbookv5/ACS_longdhbcover.html 
^See http:/ /www.stsci.edu/instruments/wfpc2/Wfpc2_dhb/WFPC2 Jongdhbcover.html 



^See htt p : / / stsdas . stsci . edu / py drizzle 
^See http:/ /stsdas. stsci.edu/pydrizzle/multidrizzle/ 
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R = 18000 at 5000A. The data were reduced using the Image Reduction and Analysis 
Facihty (IRAfJ^ software. 

The observations were carried out over several epochs from 10 days (2006 February 22) 
to 168 days (2006 July 30) after outburst. The log of spectroscopic observations is shown in 
Table [H For direct comparison with the HST imaging we focus on the spectral wavelength 
region around the [O ill] A5007 emission line, which falls in the HST F502N narrowband 
filter. 



RESULTS OF OBSERVATIONS 



3.1. Hubble Space Telescope Images 

The first epoch observations clearly show extended structure in the [O ill] A5007 line 
in the deco nvolved images, as well as in the raw images at both epochs (Figured]). As 



reported in IPaper ll . a striking feature is the double ring structure (top image in Figure [T]) 
with the major axis lying east-west with a total (peak-to-peak) extent of 360 ± 30 mas, 
in the plane of the sky, corresponding to an expansion rate of 1.2 ± 0.1 mas day~^ (and 
equivalent to T4 = 3200 ± 300 km s~^ for d = 1.6 kpc, where Vf^ is the transverse velocity) . 



Schaefer 



As an aside, we note that a much larger distance of 4.3 kpc has been suggested by 
( 120091 ). If correct, that would result in Vj^ = 8600 km s~^ which is of course very much a 



lower limit to the true expansion velocity of material and which appears to be unfeasibly 
large for a n ova outburst . We a lso note that the most recent compilation of all distance 
estimates by 



Barrv et al. 



(j2008l ) gives d = 1.4^q;2 kpc, in line with the value we assume 



^IRAF is distributed by the National Optical Astronomical Observatories, which is oper- 
ated by the Associated Universities for Research in Astronomy, Inc., under contract to the 
National Science Foundation. 
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here. More extended emission is detected above background in the deconvolved image 
with an expansion rate of 1.7 ± 0.2 m as day^^ and, when compared with values seen in 
the radio, it was suggested in 



Paper 1 



that the east-west bipolar emission seen here and in 



the radio arises from the same regions of the remnant, if the e xpansion velocities in the 



east-west direction are roughly constant after outburst (see also 



O'Brien et al. 



20081) 



The north-south (peak-to-peak) extent is 150 ± 25 mas corresponding to an expansion 
rate of 0.48 ± 0.08 n aas day~^. Results of the expansion in the early radio ob s ervati ons 



(jO'Brien et al 



2006 



Rupen et al 



20081 ) and more recently from iO'Brien et al. 



fl2008h . led 



to a north-south expansion rate of 0.77 ± 0.04 mas day ^ derived from observations over 



the first 107 days, suggesting deceleration in this direction (jPaper j ). 



In the second epoch images (bottom images on Figured]), structure is still visible 
in the [O ill] A5007 line. It has an east-west extent of 1100 ± 100 mas (peak-to-peak) 
corresponding to an expansion rate from the center of 1.2 ± 0.1 mas day~^ suggesting 
no deceleration had taken place between epochs in terms of greatest east-west extent. 
The north-south extent is much harder to determine but is approximately 460 ± 46 mas 
corresponding to an expansion rate of 0.51 ± 0.03 mas day^^. This result would imply that 
no deceleration occurred in the north-south direction between the two epochs. 

It is interesting to compare the derived extent of the remnant at each epoch with the 
expected dimensions of the red giant wind. Assuming a wind that has been expanding 
for 21 yrs (the pe riod between the last two outbursts) and a wind velocity of 40 km s~^ 



(IWallerstei 



19581 ). this implies a maximum wind radius of 2.6x10 cm. For a distance 



d = 1.6 kpc, the first epoch angular size suggests the outburst remnant has maximum 
east-west extent from the center of 4.3 ± 0.4x10^^ cm and at the second epoch the remnant 
has maximum extent 1.3 ± 0.1x10^^ cm. Therefore, by the time of the i/ST observations 
the remnant appears larger in the east-west direction than the expected size of the red giant 
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wind. 



3.2. Spectroscopy 



In Figures [2}|1] we present the spectral evolution of RS Oph as determined from low 
resolution spectra obtained f rom day 31 to day 83 after outburst. The principal emission 



lines were identified following iRosind (119871 ). Figure [2] shows the observed spectrum 31 days 
after outburst. The visual magnitude had declined to V ~ 9 by this time. The spectrum 
shows the presence of broad emission lines of H/3, H7, B.6, He I (A4471, 4713, 4922, 5016, 
5048 and 541 1), Fe ll (multiplets 27, 28, 42, 48, 49), Nlll A4640 and [Sill] A5041. What 



Rosind (119871 ) identified as the [Arx] A5535 coronal line is also present. Figure [3] shows 
spectra taken on days 64 and 65 after outburst. The magnitude had now declined to V ~ 
10. The spectra show an increased degree of excitation, with the presence of moderately 
strong coronal lines of [Arx] A5535 plus [Fexiv] A5303, and [Fex] A6374. Also present at 
this time are the nebular lines of [O ill] (A4363, 4959 and 5007) and [N ll] A5755. Figure H 
shows spectra on day 107 after outburst (V ~ 11.5, i.e. approaching minimum light). Here 
the [Fexiv] A5303 and [Arx] A5535 lines have become significantly weaker. In contrast, 
the spectrum also shows an increase in the relative strength of [O III] , and [N ll] , and [O l] 



( Rosino 


1987; 


Rosino & lijima 


1987) 



19871 ). We also note, as an asi de to ou r main discussion 



Rosind (119871 ) and others as 



that in the day 64 spectrum for example, the lines identified by 
coronal lines are blueshifted with respect to the Balmer and other lower excitation lines by 
~ 100 km s~^. This may be related to the relative principal regions of origin of these lines 
within the velocity field of the remnant. 

We used the higher resolution spectra in the region of the [O ill] A5007 emission 
line to distinguish features arising from different parts of the remnant (Figure [5]). For 
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example, during the months of 2006 February and March (bottom two panels in Figure 
El approximately ~ 10 to ~ 30 days after outburst), we identified emission lines of [O ill] 
A4959 and A5007, He I A5016 and A5048, Fe ll A5018 and [Sill] A5041. We can immediately 
see that He I A5048 and [Sill] A5041 lines are not strongly present in the spectra in 2006 
June and July (top two panels in Figure [5]). Here we interpret narrow emission lines as 
arising from the ionized red giant wind ahead of the forward shock, whereas the broad 



featu res are associated with the shocked ejecta or the shocked wind (see e.g. 



19961). 



Shore et al 



With the knowledge of the likely origin of various features in the [O ill] A5007 spectral 
region, we can start modeling the emission in this line from the expanding remnant alone 
(i.e. by ignoring emission from other elements or from the ionized wind - see Section |4] 
below). The line profile we model in detail is that of 2006 July 30 {t = 168 days) since it 
was obtained closest in time to our HST imaging (top panel Figure [5]). 



MODELING THE REMNANT 



Numerical simulations by iLloyd et al.l (1l993l ) suggested that bipolarity can be achieved 
by either (i) the outburst being spherically symmetric and the bipolarity resulting from 
subsequent interactions with an anisotropic red giant wind or {ii) the outburst itself being 
intrinsically bipolar, with material being ejected pref erentially in the po lar directions and 



Lloyd et al. 



(119931 ) found the latter 



then interacting with an isotropic ambient medium, 
to be the best fit to (the admittedly sparse) observations of bipolar structure from VLBI 
observations of the 1985 outburst of RS Oph. 



More recently. 



Sokoloski et al. 



(120081 ) suggested that the observed bipolar structure in 



the radio was not due to an intrinsically asymmetric explosion or to shaping of the ejecta by 
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circumbinary material 



disk. Other authors (ISkopal et aL 



but due to coUiraation at very early tim es, possibly by the accretion 



2008 



Banerjee et al. 



2009f ) have modeled the RS Oph 



outburst by means of Gaussian fitting to line profiles and find several components that are 
interpreted as arising from different regions of the expanding remnant; namely the central 
broad component is associated with more slowly expanding material from the waist while 
the outliers of the lines are associated by these authors with high velocity ejecta. 



Here we have used Shape fISteffen &: L6peal2006l ) to analyse and disentangle the 3-D 
geometry and kinematic structure of the expanding nebular emission in RS Oph. Shape 
was originally developed to model the complex structures of Planetary Nebulae and is 
based on computationally efficient mathematical representations of the visual world which 
allow for the construction of objects placed at any orientation in a cubic volume. It allows 
modeling of the structure and kinematics of an object to compare with observed images 
and spectra and includes parameters such as location and width of the spectrograph slit, 
seeing values and spectral resolutions. We note here that the nebular remnant was of course 
much smaller than the width of the spectrograph slit at the time of the contemporaneous 
HST observations on day 155. 

We used the information described in the previous section regarding derived velocities 
and structure of RS Oph's extended emission as initial parameters in Shape. Particles 
within Shape, each having uniform emissivity, are used to provide a simulated density and 
velocity field and can be either distributed uniformly on the surface or within the volume 
of the appropriate structure. Since the emission line of interest, [O ill] A5007, is very 
broad and extends beyond the blue cut-off of the i/ST filter, particularly for the ACS/HRC 
camera, we also applied a Doppler filter to replicate the finite extent of the HST narrow 
band filters (in the case of F502N of the ACS, the central wavelength of the filter, width 
57A is 5022A so the sensitivity to material outside the range ~ —800 to + 2600 km s~^ is 



-14- 



greatly reduced - see below). 



This modeling also allows us to determine nebular remnant inclination and relate this 
to that of the central binary (where an inclination i = 90° corresponds to the major axis 
of the bipolar structure or the orbital axis of the binary lying in the plane of the sky). 



Different autho rs have suggested that the central binary syste m has an inc. 



30° < i < 40° flDobrzvcka fc Kenvon 



1994h or 49° < z < 52° flBrandi et al. 



ination between 



20091) 



5. RESULTS AND DISCUSSION 



5.1. General Structure 



Figure [T] suggests th at the evolving resolv ed struc ture in R S Oph is bipolar and it has 



been modeled as such by 



O'Brien et al. 



(120061 ) and in 



Assuming a bipolar model, 



we then varied the inclination of the system. We also adopted a linear velocity field given 
by 



3200 r 



exp 



km s ^, 



smi ro 

where r is the distance of a particle from the center of the remnant and ro is the true 
semi-major axis. The waist was modeled so that the north-south extent was constrained to 
3.60 X 10^^ cm in line with that derived in Section [XT] above. In Shape first we introduced 
the dumbbell structure with 3 x 10^ particles distributed on the surface (as the image 
suggested the outer lobes were not centrally filled) and ran the model several times to 
produce a well sampled model spectrum. The model image was then pixelated to replicate 
the HSTACS/RRC CCD pixel size of 0.025 arcseconds. 

This simple model does not reproduce the observed spectrum. In particular, it gives 
too broad a line profile, suggesting that lower velocity material must dominate the emission 
at t = 155 days. We therefore introduced an over-density towards the center containing 
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lower velocity material. This was accomplished via an inner hour glass structure with 
semi-major axis 1.26 x 10^^ / sin i cm (for d = l.Q kpc) and with the particles distributed 
within the volume as the central structure showed no evidence of limb brightening in the 
HST image. It should be noted that the precise geometry of the inner structure may not be 
as well constrained as that of the outer dumbbell. However, we have for example explored 
a simple spherical geometry for this inner region, but the fits to the observed spectra were 
then much poorer in terms of the gross line profile (inverted U-shape) and finer detail (no 
double central peak). It should also be noted that this inner structure is smaller than the 
estimated extent of the red giant wind at this epoch. 

The model system thus contains two structures with the central hour glass having four 
times more emitting particles than the outer dumbbell structure (Figure [6]); this ratio is in 
approximate agreement with the estimated ratio of flux from the inner region of the image 
to that in the lobes (factor of 4; likely to be a lower limit for this ratio as some of the flux 
in the central region will be contributed by the innermost parts of the outer structure). 

The existence of a distinctly two component structure is compatible with early 



infrared interferometry (jChesneau et al. 



2007 



Chesneau 



20081 ). and r adio observations o f 



a central peak (thermally dominated) and outer lobes (non-thermal; iTaylor et al. 



Sokoloski et al. 



1989 



2008 



Eyres et al. 



20091 ). In addition, Vaytet et al. (in preparation) have 
modeled the early time X-ray spectra of RS Oph and found that they also required a 
density enhancement towards the central regions to give a high enough absorbing column. 



5.2. System Orientation 



Figure [7] shows the results of modeling using the outer dumbbell and inner hour glass 
structures described above. The model image reproduces well the general morphology seen 
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in the HST image, including that of the west side appearing less prominent than the east. 
This asymmetry is an observational effect, due to the finite size of the HST filter F502N 
(as illustrated in Figure [3 see also Section 4 above) and immediately tells us that the 
western lobe is approaching the observer. Indeed, any apparent disagreements between the 
model and observed images can largely be attributed to artefacts of the PSF subtraction 
procedure. We note that the spectrum shows excess emission on the blueward side (top 
Figure [5]). This could be attributed to Fe ll A5001.91, although this line is not observed in 
spectra at earlier times (Figure [5]). 



Although radio observations led 



O'Brien et al 



(120061 ) to model the RS Oph remnant 



as a bipolar structure, the deduced orientation of the system differs from that found from 



our results in the optical. 



O'Brien et al 



(I2OO6I ) modeled the orientation with the east lobe 



closest to the observer since that is where a second component appears first in their radio 
observations at early times. However, our observations suggest that the west lobe must be 
approaching the observer. This is in agreement with results from infrared interferometry at 
early times and it is suggested that the early emergence of the eastern-most radio lobe may 



be due to obscuration 



(IChesneau et al. 



2007 



jy a flared disk of circumbinary material lying in the orbital plane 



Chesneau 



20081 ). Even at the epoch of the first VLBI observation 
at f = 13.8 days, this material has to be external to the binary and would therefore arise 
from an anisotropic red giant wind. We also find that the system has a position angle on 
the sky of 85 ° in lin e with that found from the radio observations of the 1985 outburst by 



Porcas et al. 



mm . 



5.3. System Inclination 

Shape also allows us to determine the inclination of the system to the line of sight. 
If the nebula's orientation is linked to that of the central binary, and with the knowledge 
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that RS Oph is not an eclipsing system, we calculated the upper limit of i to be used in 
our exploration of the incli nation angle. For an e clipsing binary with an M2 III giant 



with a radius of 67^11 R© (IDumm &: Schild 



19981 ) and a separation between the giant 



and WD of ~ 1.5 AU (if for example I 


^wn - 


f Mnn = 2Mm anc 


Porb 

)) we 


Dobrzvcka & Kenvon 


1994; 


Fekel et al. 


2000; 


Brandi et al. 


2009 



occur if i > 78±3 degrees. A lower limit for the range over which we explored the model 
was constrained by velocities being far too high to arise from a nova explosion; e.g. at i = 
20°, Vexp = 9300 km s~^. We therefore ran the model for 20° < i < 81° and determined 
each model's respective fit by comparing the observed and model spectra. 

We find the best value for the inclination of the remnant to be i = 39tio degrees (Figure 
[HI the errors were determined in the standard way by taking the AX'^ and calculating 
th e la level). This is in g ood a greement with the inclination for the central binary found 



by 



Dobrzvcka &: KenyonI (119941 ) o 



determination of 49° < i < 52° by 



30° < i < 40°, th ough less so with the more recent 



Brandi et al 



20091 ). For illustration. Figure [7] also shows 



the synthetic spectra derived at the error limits quoted for the inclination where clearly the 
synthetic line profiles do not reproduce as well the observed profiles as does the best fit 
value of i = 39°. We note that lines of Fe ll A5001.91, He I A5015.68 and Fe ll A5018.44 have 
been removed from the observed spectrum prior to model fitting, following the procedure 
outlined in Section [3|, to ensure that our spectral fits are only to the [O ill] A5007 line itself. 



5.4. Second Epoch 

We also modeled the second epoch HST observation at t = 449 days. Using the built-in 
time evolution option in Shape, we first explored a simple linear expansion from the first to 
the second epoch using the results found in Section 4. In this case, the effective wavelength 
(5OI2A) being much less offset from the rest wavelength of the [O ill] A5007 line itself for 
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the same filter in tlie ACS means tliat tlie asymmetry in the resulting image introduced 
by the Doppler cuts is not as marked as it was for the first epoch (ACS) image. Figure M 
shows the results. 

Little resemblance is evident between the observed and modeled images in this case. 
This is an important result because it tells us that the second epoch is not just a linear 
expansion of the first, at least in the central regions. We then allowed the dumbbell 
structure to expand linearly but the inner hour glass was kept the same size as in the 
first epoch (right hand image in Figure [9]). We may understand the linear expansion of 
the outer structure if, as we calculated above, it had already cleared the pre-existing red 
giant wind by the first epoch. Alternatively, the apparent constant expans i on rat e may 



Sokoloski et al. 



(120081 ). The 



be consistent with powering by narrow jets as proposed by 
apparent constancy of the inferred size of the innermost hour glass structure is more 
difficult to understand. However, we should be aware of over-interpretation of the poorer 
image from the second epoch, plus the lack of contemporaneous optical spectroscopy to 
aid our modeling. Although the model reproduces some other features of the second epoch 
observations, for example the gap in the east lobe due to the applied Doppler filter, further 
investigation is required. 



6. SUMMARY 

We have reported combined HST imaging, ground-based spectroscopy and detailed 
modeling of the RS Oph nebular remnant following the 2006 outburst. We find that: 

1. The remnant of RS Oph can be understood as a bipolar structure with two co-aligned 
components. In order to satisfactorily model both the image and spectra from day 
155, we require the presence of a high density (low expansion velocity) central region 
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and a more extended (higher expansion velocity) less dense region in the form of two 
lobes. This morphology is in agreement with that inferred from observations in the 
radio and at X-ray wavelengths, plus infrared interferometry at early times. 

2. The asymmetry observed in the first epoch HST image is due to the finite width of 
the HST filter. This in turn implies that the west lobe is approaching the observer, in 
agreement with the results of infrared interferometry at early epochs. 

3. We are able to determine the inclination of the remnant as 39 lin degrees. This is 
compa rable to estimates of the inclination of the central binary of 



Dobrzvcka fc Kenyon 



( 1l994j ). but lies outside the range of i from recent estimates by iBrandi et al 



mm 



4. Modeling of the second epoch HST image implies that the outer dumbbell structure 
underwent linear expansion; however, there is more evidence of deceleration for the 
central (hour glass) region. 
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Fig. 1. — HST images of RS Oph through the F502N filter. Top left - Raw image from the 
first epoch observations (t = 155 days after outburst) using the ACS/HRC. Top right - PSF- 
subtracted and deconvolved image using a TinyTim PSF, showing a double ring structure. 
Bottom left - Raw image from the second epoch observations (t = 449 days after outburst) 
using WFPC2. Bottom right - PSF-subtracted and deconvolved image using a TinyTim 
PSF, again showing evidence of a double ring structure (note that the "jet-like" feature on 
the WFPC2 deconvolved image is an artifact caused by bleeding on the CCD chip). North 
is up and east is to the left in all images. 
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Fig. 2. — Spectral observations of RS Oph on 2006 March 15 (t = 31 days after outburst - 
see Table [1]) . 
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Fig. 3. — As Figure [21, but for 2006 April 17 and 18 (t = 64 — 65 days after outburst). Note 
that the day 64 spectrum has been offset in flux by +1 dex for clarity. 
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Fig. 4.— As Figure El but for 2006 May 29 (t = 107 days after outburst). 
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Fig. 5. — RS Oph [O ill] A5007 spectral evolution at different epochs after the 2006 outburst. 
From bottom to top, 2006 February, March, June and July (dates as shown). Here the 
evolution of [O ill] A5007 can be tracked first as arising from the ionised wind of the red 
giant (narrow lines in February and March) and then at later times (June and July) it is 
more likely related to shocked ejecta and/or shocked wind (broad lines). See Table [1] for 
more details of spectroscopic observations. 
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Fig. 6. — Dumbbell (outer) and hour glass (inner) structure that is used to replicate the RS 
Oph remnant. 
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Fig. 7. — Top - Synthetic image without the Doppler filter apphed (left), enlarged HST/ ACS 
image at t = 155 days (middle) and synthetic image with Doppler filter applied (right). The 
synthetic images are from a combined outer dumbbell and inner hour glass structure inclined 
at 39 degrees to the line of sight and with a position angle of 85°. Bottom - Best fit synthetic 
spectrum (black; left) is overlaid with the observed spectrum (red; the emission lines of Fe ll 
A5001.91, He I A5015.68 and Fe ll A5018.44 associated with the ionized wind ahead of the 
forward shock have been removed). Model spectra resulting from inclination angles at the 
la error limits are also shown (right) where the dotted line is for i = 40° and the dashed 
line is for i = 29°. 
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Fig. 8. — AX^ as a function of inclination angle. This reveals that the best-fit value for 
the inclination of the system is 391]^q degrees {la errors; 3a errors are +8° and -14°). The 
inclination was varied in 1° steps and then, when close to the minimum, in 0.1° steps (note 
that we explored inclination angles up to 81° as suggested in section 5.3, but the models 
continued to provide a very poor fit to the data at such high inclinations). 
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o:'4 

Fig. 9. — Left - results from a linear expansion of the first epoch model. Middle - HST 
WFPC2 image t — 449 days. Right - result from a model in which the outer dumbbell 
structure was expanded linearly but the central hour glass was kept the same size as in the 
first epoch. The model images have a Doppler cuts applied at around —1600 km s~^ and 
+2300 km s~^ as appropriate for the WFPC2 filter. North is up and east is to the left. 
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Table 1. Log of optical spectral observations. 



Date (days 
after outburst) 


Instrument 


Total Exposure 
time (s) 


R 


zUUD reoruary zz [iO) 


Echelle 




iDUUU 


zUUd reoruary zd (iij 


Echelle 


4zo 


iDUUU 


zUUo reoruary z4 


Echelle 




ioUUU 


2006 February 25 (13) 


Echelle 


435 


16000 


2006 March 14 (30) 


Echelle 


600 


16000 


2006 March 15 (31) 


Boiler & Chivens 


90 


2355.2 


2006 April IT'^ (64) 


Boiler & Chivens 


150 


1668.6 


2006 April 18* (65) 


Boiler & Chivens 


140 


3149.6 


2006 May 29 (83) 


Boiler & Chivens 


40 


1525.9 


2006 June 18 (126) 


Echelle 


1830 


16000 


2006 June 19 (127) 


Echelle 


990 


16000 


2006 July 30 (168) 


Echelle 


420 


16000 



*Observations taken at the Observatorio Astrofisico Guillermo Haro, 
otherwise from Observatorio Astronomico National en San Pedro 
Martir. 



